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Abstract: Deregulated immune response to gut microflora in genetically predisposed 
individuals is typical for inflammatory bowel diseases. It is reasonable to assume that 
genetic association with the disease will be more pronounced in subjects with early onset 
than adult onset. The nucleotide-binding oligomerization domain containing-2 gene, 
commonly involved in multifactorial risk of Crohn's disease, and interleukin 10 receptor 
genes, associated with rare forms of early onset inflammatory bowel diseases, were 
sequenced in an early onset patient. We identified a novel variant in the N0D2 gene 
(C.2857A > G p.K953E) and two already described missense variants in the ILIORA gene 
(S159G and G351R). The new N0D2 missense variant was examined in silico with two 
online bioinformatics tools to predict the potentially deleterious effects of the mutation. 
Although cumulative effect of these variations in the early onset of the disease can be only 
hypothesized, we demonstrated that family information and in silico studies can be used to 
predict association with the disease. 
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1. Introduction 

Inflammatory bowel diseases (IBDs) are chronic inflammatory disorders of intestine that result 
from a deregulated mucosal immune response in genetically predisposed individuals [1]. It is a 
widespread pathological condition with increased incidence in Western populations, affecting 
approximately one out of 250 people. The two major forms of IBD, Crohn's disease (CD) and 
ulcerative colitis (UC), are usually quite distinct. Although the age of onset of IBD is usually in early 
adulthood, the pathology could also occur in childhood during the first years of life. Early onset IBD 
patients often present overlapping features of CD and UC, making a definite diagnosis difficult. 
Various epidemiological studies suggest that patients with an early onset diagnosis usually manifest a 
more serious and aggressive course, tending to a very fast development [2]. Evidence supports the 
hypothesis that genetic predisposition plays a more substantial role in pediatric patients, who manifest 
a more extensive and severe course of the disease than subjects with adult onset. 

Several genes, identified by several genome-wide association studies (GWAS), have been 
associated with IBD predisposition, both in adult and early onset. In addition, a number of monogenic 
immune disorders have been found in a small, but increasing, portion of early onset IBD, suggesting 
that common genetic variants and rare monogenic defect can play a synergic and crucial role in the 
development of the disease. 

Considering that GWAS could only explain about half of the genetic heritability of IBD, it has been 
suggested that the additive role of multiple variants and/or rare "private" genes not identified by 
GWAS may account for the so-called "missing heritability"; this could be even more likely for early 
onset IBD, in which higher genetic risk is expected. 

The patient we are going to describe is a child (present age: seven years) presenting a very early 
onset ulcerative colitis (VEO-UC), cared for at the Pediatric Gastroenterological Service of the 
Institute for Maternal and Child Health IRCCS "Burlo Garofolo" (Trieste, Italy): at 18 months, 
pancolitis with mucosal and intracryptic localization was diagnosed. No extra-intestinal manifestations 
and/or associated diseases were present. At colonoscopy, rectal inflammation with hemorrhagic ulcers 
and mucosal edema were evidenced. The colonic mucosa was very fragile and tended to bleed easily. 
Intestinal biopsies showed chronic inflammation with lymphocytic aggregates. The search for 
granuloma and viral inclusions results as being negative. The disease showed a severe course in spite 
of the use of glucocorticoids, azathioprine and thalidomide administration. A good clinical response 
was achieved only with a biological therapy with infliximab (anti-TNF-a antibodies). In order to assess 
monogenic causes of early onset inflammatory colitis in this patient, we analyzed both subunits alpha 
and beta of the interleukin-10 receptor (ILIORA and ILIORB), as well as nucleotide-binding 
oUgomerization domain containing 2 (N0D2), since these genes are known to be associated with a 
higher risk for CD. 

2. Results and Discussion 

2.1. Results 

We found 18 variants in our patient, five in the N0D2, four in the ILIORA and nine in the ILIORB 
genes. All variants localized respectively at the 5' and/or 3' untranslated, intronic and coding regions 
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(Table 1). Among the variants identified in N0D2, four are known variants, and one, is a novel 
missense variant at the exon 9 (c.2857A > G p.K953E) present in heterozygosis (Figure IB). Within 
the three variants in the coding sequence of ILIORA, two missense variants, both present in 
heterozygosis, rs3135932 (c.475A > G p. S159G) and rs2229113 (c.l051 G > A p.G351R), have 
already been described in the literature. 

Table 1. Genetic variant in the analyzed genes, relative genotyping and localization. 
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The three-missense variants were searched for in both parents to check their pattern of inheritance. 
The mother carried the three variants (K953E, S159G and G351R) observed in the patient, while the 
father results heterozygous only for the G351R variant (Figure 1). The new N0D2 missense variant 
was searched for in 60 anonymous healthy Italian donors to confirm if it was not a polymorphism, but 
a true mutation: no variant has been found. 

All variants were already described in the common SNPs Database (Single Nucleotide Polymorphisms 
at http://www.ncbi.nlm.nih.gov/snp/), as polymorphic in the general population, except for the new 
N0D2 variant, for which we performed an in silico analysis with PolyPhen2 (Polymorphism 
Phenotyping v.2) and SIFT (Sorting Intolerant from Tolerant). Both online tools showed a close 
agreement in the scores, predicting this new missense variant, K953E, as probably damaging. Several 
alignments in different species showed that this is a conserved region (Figure 2). 
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Figure 1. Sequence chromatograms from patient and parents of the ILIORA and N0D2 
coding regions. (A) In the ILIORA exon 4 and exon 9, the arrows indicate the nucleotide 
substitution, C.475A > G and c. 1051 A > G, consisting, respectively, in the amino acid 
substitutions, S159G (A/G heterozygous patient and mother, A/A wild type father) and 
R351G; (B) in the N0D2 exon 9 sequence, the c.2857 A > G substitution consisted in an 
amino acid substitution, K953E (A/G heterozygous patient and mother, A/A wild-type father). 
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2.2. Discussion 

In our case report study, we analyzed a set of three genes, described as strongly associated with 
IBD [3,4], in a child affected by severe early onset inflammatory colitis. Among the 18 variants, three 
missense variations, the new N0D2 mutation (K953E) and the known ILIORA missense variants 
(S159G and G351R), caught our attention: we hypothesized that these three variants might together 
contribute to an increased risk of developing early onset IBD [5]. Variants in combined heterozygosis 
on different genes may have a cumulative effect and contribute to the disease phenotype, probably due 
to the dysregulation of different pathways. Digenic contribution to multigenic diseases is poorly 
investigated by GWAS studies. However, recent data showed that digenic or multigenic inheritance 
might be underrated mechanisms explaining familial recurrence of multifactorial diseases [6]. 

ILIORA (llq23) and ILIORB (21q22) are attractive inflammatory bowel disease candidate genes: 
they encode for the alpha and beta chain, respectively, that form a tetrameric cell bound structure able 
to bind interleukin (IL-10). The IL-10 pathway is very important in regulating the intestinal 
inflammation, and mutations in the receptor cancel the IL-10 immunomodulatory signal, which is 
strongly associated with an intestinal hyper-inflammation [7]. Although there are no data about the 
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association between these variants (S159G and R351G) and the early onset of the disease, several 
studies have shown a loss of function due to these variants [8,9]. 

N0D2, localized at chromosome 16q21, codify for a protein that belongs to the family of 
intracellular NLR (NOD-like receptors), able to recognize microbial components and to stimulate an 
inflammatory response through the activation of NF-kB. Moreover, three N0D2 mutations 
(rs2066844, rs2066845 and rs5743293) represent the main genetic factor causing susceptibility to CD. 
These mutations occur in the LRR domain or in the neighboring domain, suggesting that they alter the 
receptor capacity to recognize the bacterial component [10]. 

Although the new K953E mutation in N0D2 has not been previously reported, it is likely to disrupt 
N0D2 protein function for several reasons: (1) the mutation is a disfavored substitution of an evolutionarily 
conserved leucine residue repeat; (2) the mutation occurs in the LRR domain, suggesting that it could 
alter the recognition of the bacterial agent; (3) the amino acid, K, is highly conserved and the 
resulting amino acid change was predicted to be deleterious by both SIFT [11] and PolyPhen2 
software [12,13] (Figure 2). However, also the mother of the IBD affected child was heterozygous for 
this mutation, not presenting, nevertheless, any disease signs. Therefore, even if in silico analysis 
considered the K953E N0D2 mutation as harmful, we are aware that further studies are needed to 
elucidate the importance of this variant in the setting of early severe colitis. 

We can speculate that the new mutation in the N0D2 gene described in our early-onset UC patient 
could alter the normal function of the receptor that plays a central role in the pro-inflammatory 
response and bacterial sensing. The other two variants in the subunit alpha of the ILIOR are responsible 
for a decreased activity of the receptor. Therefore, our assumption is that, although these variations are 
located on two different genes and apparently belong to different molecular pathways, they could 
somehow determine a cumulative dysfunctional effect, causing themselves the peculiarities of this 
form of ulcerative colitis. 

We are aware that these are just hypotheses that should be confirmed by functional studies, lacking 
in this case report and representing a limit of our study. Moreover, the presence of N0D2 (K953E) and 
ILIORA variants in the healthy mother suggests that the associated risk may not be so high. This brings 
us back to the concept of multifactorial disease in which also the non-genetic factor should be 
considered: for this reason, we can therefore hypothesize that the disease exhibited by the child could 
be given by the combination of the three variants with an additional factor, such as diet, antibiotic 
consumption and microbiome, and other not yet investigated genes. Thus, further analyses, including 
functional studies to verify the true effect of N0D2 mutation and ILIORA, should be performed, 
recruiting a larger number of patients in order to better understand the genotype -phenotype correlation, 
in relation to an additional factor or an additional gene. 
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Figure 2. Bioinformatics analysis results. (A) Multiple alignment of the amino acid 
sequence of N0D2 protein in seven species showed that this is a conserved region; 
(B) PolyPhen2 (Polymorphism Phenotyping v. 2) analysis predicting the probably 
damaging impact of the K953E substitution with a score of 0.999. 
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3. Experimental Section 

3.1. DNA Extraction and Genotyping 



Genomic DNA was extracted from 1-2 mL EDTA-anticoagulated blood using an EZl DNA Blood 
Kit (QIAGEN, Milano, Italy), according to manufacturer's instructions. PGR amplification was 
performed for the entire coding sequence, flanking regions and for 5' and 3' untranslated region of the 
N0D2 and III OR genes, using the KAPA 2G Fast Hot Start Readymix (RESNOVA, Roma, Italy) and 
using a 2720 Thermal Gycler (Applied Biosystems, Foster Gity, GA, USA). Amplification products 
were directly sequenced by the Sanger method using an ABI PRISM 3130XL automated DNA 
sequencer (Applied Biosystems, Foster Gity, GA, USA). Sequences were analyzed with Seqman II 
Software (DNASTAR I Lasergene, 7.0, Madison, WI, USA). While primers used for amplification of 
ILIO receptor genes (ILIORA and ILIORB) were designed according to the indications of 
Glocker et al. [3]; N0D2 primers used both for PGR and sequencing were designed using Primer Blast 
based on the relative sequence deposited in GenBank: NM_022162 (NOD2/CARD15). Each variant 
has been identified and verified in databases by performing a nucleotide BLAST (Basic Local 
Alignment Search Tool) server (http://blast.ncbi.nlm.nih.gov/). 

3.2. In Silico Analysis 

The bioinformatics analysis of missense variants was performed with two online tools that can 
predict the possible impact of an amino acid substitution on both the structure and function of a human 
protein using straightforward physical and comparative considerations: PolyPhen2 (version 2.2.2, 
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PMID: 20354512 and 23315928, Boston, MA, USA), (http://genetics.bwh.harvard.edu/pph2/) and SIFT 
(version 1.03, PMID: 11337480 Seattle, WA, USA), (http://sift.jcvi.org). 

4. Conclusions 

Finally, in the context of this kind of genetic disease, we would like to emphasize that genetic 
analyses of many genes can open new scenarios that could be unraveled only by large studies and 
robust biologic and bioinformatics evaluations. For this reason, a new generation sequencing (NGS) 
approach must be taken into account to seek an alternative genetic cause that may contribute to the 
disease phenotype. 
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